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g ■ Abstract 

The anomalous weak dipole moments of the r lepton are measured in a data sample collected 

■ by ALEPH from 1990 to 1995 corresponding to an integrated luminosity of 155 pb _1 . Tau leptons 

^ • produced in the reaction e + e~ — > t + t~ at energies close to the Z mass are studied using their 

Q-i! semileptonic decays to tt, p, ai — > tt2tt° or a\ — > 3it. The real and imaginary components of both 

c-| ! the anomalous weak magnetic dipole moment and the CP-violating anomalous weak electric dipole 

; moment, Re/i T , im/i r , Re<i T and lmd T , are measured simultaneously by means of a likelihood fit 

k> ' built from the full differential cross section. No evidence of new physics is found. The following 

^2 : bounds are obtained (95% CL): \Re/i T \ < 1.14 xlO" 3 , \Im/i T \ < 2.65 x 10~ 3 , |Rerf T | < 0.91 x 10~ 3 , 

. and |Imd T | < 2.01 x 10~ 3 . 
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1 Introduction 



The anomalous weak dipole moments of the r lepton are the tensorial couplings of the Zt + t~ 
vertex. They are zero to first order in the Standard Model (SM). Two types of anomalous weak 
dipole moments can be distinguished: the magnetic term \i T and the CP-violating electric term 
d T . Here, both the real and the imaginary components of each anomalous weak dipole moment 
are explored, i.e. Re/z r , lm.^ T) Ked T and Im<i r . Radiative corrections in the SM provide nonzero 
predictions for fj, T and d T [1, 2] which are below the present experimental sensitivity. This opens 
the possibility to look for deviations from the SM. 

There have been many searches for the CP-violating anomalous weak electric dipole moment 
of the r since the beginning of LEP [3-5]. In addition, limits on the anomalous weak magnetic 
dipole moment were obtained more recently [5]. 

In this analysis, the previous ALEPH result on Reci r [3] is updated, and Re/i T , Im/i r and 
Im d T are determined for the first time in ALEPH. The data sample was collected with the ALEPH 
detector from 1990 to 1995 at energies around the Z resonance and corresponds to an integrated 
luminosity of 155 pb _1 . Tau leptons are generated in the reaction e + e~ — > t + t~ at LEP. The 
method to extract the anomalous weak dipole moments is based on a maximum likelihood fit to 
the data taking into account all r spin terms explicitly, including correlations. This is the first 
time that the complete differential cross section for the production and decay of the r leptons is 
considered to estimate the r anomalous weak dipole moments. The most important semileptonic 
decays are used: n, p, a\ — ■> 7r27r° and a\ — > 3n. The r spin information is recovered using optimal 
polarimeters which are different for each decay. The selection and particle identification make use 
of tools already developed in previous analyses [6-9] . 

The text is organized as follows. The theoretical framework is introduced in Section 2. The 
most important ALEPH subdetectors for this analysis are covered in Section 3. The data analysis 
procedure is explained in Section 4, emphasizing the new features of the analysis. The more 
relevant systematic uncertainties are then discussed in Section 5. The results and conclusions are 
presented in Section 6. 



2 Theoretical framework 



2.1 Production cross section 



The currents assumed for photon and Z exchange in e + e — > t + t production are 



-<m(7) 



iQfej^, with / = e, r , 



T^ (z) = ie[v e -f-a e -f 



jVt(Z) 



75] 



2m T 



2m T 



lh<r~q v 



(1) 



where Qfe is the fermion charge; a e , a T , v e and v T are the axial vector and vector couplings of the 
SM; fj, T and d T are the anomalous weak magnetic and anomalous weak electric dipole moments 
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of the r. In the previous expression both anomalous weak dipole moments are dimensionless 
quantities. However, the anomalous weak electric dipole moment is often quoted in the literature 
in units of ecm, by defining the contribution of this dipole moment to the current as id r ^^a^ v q u . 
These different notations are related by the conversion factor e/2m T = 5.552 x 10 _15 ecm. 

Using the currents in Eq. 1, the differential cross section can be expressed as [10] 

-, d<J ' (si,s 2 ) = R o+ J2 R no- s i+ RovS V 2+ R^uSiS V 2 . (2) 

« cost V /i=l,3 v=l,3 /i,f=l,3 

The R^y terms are functions of the fermion couplings and of the r production angle 9 T ; s\ and s 2 
are unit vectors chosen as the quantisation axes for the spin measurement of the r + and the r~, 
respectively, in their corresponding rest frames. 

The following reference frame has been chosen: the z axis is in the outgoing r + direction and 
the incoming e + is in the yz plane. The x component is therefore normal to the production plane. 
The y component is called transverse. 

Several R^ v terms have been already measured by ALEPH. Defining {R jXV ) ± = (R^ v ± -R„ M ), 
these terms are the following: 

- R 00 = da/ dcos9 T [11], 

- (i?o3)+/-Roo — P T (cos9 T ) is the longitudinal polarisation of the r [9], 

- R22/R00 = —R11/R00 are the transverse-transverse and normal-normal spin correlations [12], 

- (i?2i)+/-Roo are the transverse-normal spin correlations [12]. 

The R^ u terms most sensitive to Re/z T , Im/z T , Re<i T and Im<i T are presented below. These 
terms are obtained from Ref. [10] after some algebra. The SM contributions are separated from 
the anomalous (anm) contributions. 

Re ji T : 



2 1 

(-^02)4. oc — sin6> T \v T \ 2 Re{y e a* e ) H sin# T cos6* T (|a e | 2 + \v e \ 2 )Rz(y T a* T ) 

SM It It 

2 1 |„. |2\ 



(-^02)4. oc It sin6* T cos9 T (\a e \ +\v e \ )i?e(a T/ u* 



anm 

2 



+ ^(Tr+jO sin 9 T Re(v e a* e )Re (v T fx* T ) + 2 7t sin 9 T Re(v e a* e ) | // T | 2 (3) 

7r 
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(-^32)4. gM oc — sm 9 T Re(v e al)Re(v T a*) -\ sin# r cos9 T (|a e | 2 + |t> e | 2 )|tv| 2 



7 2 + 1 

(^32)4. oc — sin9 T cos9 T (\a e \ 2 + \v e \ 2 )Re(v T fj,* T ) 

"r anm -y^ 

+ 2j T sin 9 T Re{y e a* e )Re{a T jj* T ) + 7 T sin 9 T cos9 T (\a e \ 2 + \v e \ 2 )\/i T \ 2 (4) 
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Im \i T : 



oc — sin# r cos6 T (\a e \ 2 + \v e \ 2 )Im(v*a T ) 

SM ry T 

oc 7 T sin# T cos6> T (|a e | 2 + \v e \ 2 )Im{a T jj* T ) 
2(7^ - 1) ^ 



+ 



7r 



■ sin 9 T Re(v e al)Im(v T /j,*) 



(5) 



RecL 



(Roi) 

(Rm). 



sm oc — sin 9 T Re{v e a* e )Im{v*a T ) 



oc 



It 



sin# T cos# T (|a e | 2 + \ve\ 2 )Irn(v T fj,*) 



+ 2'-f T sin 6 T Re{v e a f e )Im{a T ii' T ) 



(6) 



(#oi) 
(Roi). 



SM 



oc 

oc — 7 T sin^ T cos# T (|a e | 2 + |t> e | 2 )i?e(a T <i*) 

— 27 T sin 6 T Re(v e a* e ) [Re{v T d* T ) + Re(n T d* 



(7) 



Im gL 



(^3l)_ 



SM 



oc 

oc — 7 T sin# T cos# r (|a e | 2 + |f e | 2 ) [Re(v T d* T ) + Re(fi T d*)} 
— 27 T sin# T Re(v e a* e )Re(a T d*) 



(8) 



(i2 32 ). 
(# 32 )_ 



SM 



oc 

oc 7 r sin# r cos# r (|a e | 2 + \v e \ 2 )Im(a T d* T ) 

+ 27 T sin 6 T Re{v e a* e ) [Im{y T d* T ) + Im(fi T d*)] 



(9) 



(i?02). 
(i?02)_ 



SM 



oc 

oc 7 r sin^ r cos9 T (|a e | 2 + |f e | 2 ) \Im{v T d* T ) + Jm(/i r rf*)] 
+ 27 r sin 9 T Re(v e a* e )Im(a T d*) 



(10) 



Taking into account that a; ^> vi, the terms can be ordered in sensitivity, and the most sensitive 
term for each anomalous weak dipole moment is presented first. The quantity 7 r is computed as 
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v /s/2m T . The photon exchange terms are omitted from these expressions for simplicity, although 
they are taken into account in the final results. 

The anomalous weak dipole moments are extracted including all R^ v terms in a maximum 
likelihood fit. In this analysis (i?3i) + , the most sensitive term to Im/i T , is used for the first time 
as proposed in Ref. [13]. The terms (#02)+, (-Roi)_ and (#32) _ were previously used in other 
measurements of the anomalous weak dipole moments. 



2.2 Tau decay 

For each r decay mode, the differential partial width of a polarised r is written as 

dT(s) = W{\ + h ■ s) dX , (11) 

using the expressions for W and h from the TAUOLA Monte Carlo program [14]; W is the 
differential partial width of an unpolarised r, and the h vector is the polarimeter of the particular 

— * 

decay mode considered. Both W and h depend on the four-momenta of the final state particles 
in the r rest frame, and they are different for each decay topology. The simplest expressions are 
those of the r decay into tt. In this case, h n is proportional to the tt momentum in the r rest frame 
and W n is a constant. In the above equation, X is a set of independent variables describing the 
full decay configuration. The number of elements of the set depends on the number of particles 
in the final state. The set X w denotes the set of variables expressing the tt direction in the r rest 
frame. 

— * 

In this analysis, the expressions for W and h allow the spin information for all the r decays 
to be recovered optimally. 



2.3 The full differential cross section 



Once the production cross section and the partial decays are introduced, the full differential cross 
section of e + e~ — > t + t~ — > x\x 2 v T v T is built following Refs. [15, 16], namely 



do 



dcos9 T dXidX 2 



r r 

T T 



/i=l,3 u=l,3 fi,u=l,3 



(12) 



In this equation T T is the total r width, X 1 and X 2 are the sets of independent variables, and hi 

— * 

and h 2 are the polarimeters for the decay of the r + and the r~, respectively. 

With the definitions = R^/Rq , H" = Wh»/T T (p, v = 0, . . . , 3), and h° = 1, the 
likelihood of an event with the final state topology ij is written as 

L ij (iJi T ,d T \9 T ,Wi,coa9h 1 ,<i>h 1 ,W2, cos 9 h2 ,<f) h2 ) = (13) 

R^r,d T ,9 T )Hj i (W 1 ,cos9 hl ,(f) hl )Hj(W 2 ,cos9 h2 ,(f) h2 ) . 

/u,i/=0,...,3 

The indices (i, j) refer to the decay mode of each r, with i,j = tt, p, 7r27r°, 3ir, and the quantities 
(Wi, cos^, 0hj and (W 2 , cos9h 2 , <f>h 2 ) are the observables related to the decay of the r + and 
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the r~, respectively. The angles (9^, (f>hu 9h 2 , 4>h 2 ) are the polar and azimuthal angles of the 
polarimeters of each r, in the reference frame introduced in Section 2.1. The above likelihood is 
also a function of the centre-of-mass energy. The distributions of the hemisphere observables W 
and cos 6^ are presented in Figs. 1 and 2. 

This likelihood fulfills the normalisation condition 

^2 J Lij(n T ,d T \0 T , W 1 ,cos9 hl ,(f) hl ,W 2 ,cos9 h2 ,(f) h2 )dX 1 dX2 = 1. (14) 

ij 

This is the integral over all possible decay parameters and all possible decay topologies ij for a 
given e + e _ — > t + t - event. The normalisation is such that the likelihood depends only upon the 
net spin polarisation of the produced r pairs, and not upon R 0Q = da jd cos 9 T . 



3 Apparatus 

The ALEPH detector is described in detail in [17] and its performance in [18]. 

Charged particles are measured with a high resolution silicon vertex detector (VDET), a 
cylindrical drift chamber (ITC), and a large time projection chamber (TPC). The momentum 
resolution in the axial magnetic field of 1.5 T provided by a superconducting solenoid is Ap/p 2 = 
0.6 x 10^ 3 (GeV/c)~ 1 for high momentum tracks. The impact parameter resolutions for high 
momentum tracks with hits in all three subdetectors are a ri p = 23 /im and a z = 28 /im. 

The tracking devices are surrounded by the electromagnetic calorimeter (ECAL), which is a 
highly segmented lead/proportional- wire-chamber calorimeter. The calorimeter is read out via 
cathode pads arranged in projective towers covering 0.9° x 0.9° in solid angle and summing the 
deposited energy in three sections of depth. A second readout is provided by the signals from the 
anode wires. The energy resolution is a/E = 0.009 + 0.18/^/E(GeV). 

The ECAL is inside the solenoid, which is followed by the hadron calorimeter (HCAL). 
Hadronic showers are sampled by 23 planes of streamer tubes giving a digital hit pattern and 
an analog signal on pads, which are also arranged in projective towers. This calorimeter is used 
in this analysis to discriminate between pions and muons. Outside the HCAL there are two layers 
of muon chambers providing additional information for ji identification. 



4 Data analysis 

4.1 Selection and decay classification 

Events from Z — > t + t~ are retained using a global selection in which each event is divided into 
two hemispheres along the thrust axis. The selection is that used in the ALEPH measurement 
of P T (cos9 T ) with the r direction method [9]. Additional information can be found in [8] and 
references therein. 

The charged particle identification is based on a likelihood method which assigns a set of 
probabilities to each particle. A detailed description of the method can be found in Refs. [6, 7]. 
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The probability set for each particle is obtained from (i) the specific ionisation dE/dx in the TPC, 
(ii) the longitudinal and transverse shower profiles in ECAL near the extrapolated track and (Hi) 
the energy and average shower width in HCAL, together with the number out of the last ten of 
HCAL planes that fired and the number of hits in the muon chambers. 

The photon and it reconstruction is performed with a likelihood method which first 
distinguishes between genuine and fake photons produced by hadronic interactions in ECAL or by 
electromagnetic shower fluctuations [8]. All photon pairs in each hemisphere are then assigned a 
probability of being generated by a tt°. High energy 7r° with overlapping showers are reconstructed 
through an analysis of the spatial energy deposition in the ECAL towers. All the remaining single 
photons are considered and those with a high probability of being a genuine photon are selected as 
7t° candidates. Finally, photon conversions are identified following the procedure described in [8] . 
They are added to the list of good photons and are included in the 7r° reconstruction. 

The r decay classification depends on the number of charged tracks and their identification, and 
on the number of reconstructed 7r°. It follows the classification for the measurement of P T (cos8 T ) 
with the r direction method, described in [9] and the references therein. 

The r selection efficiencies and the background fractions for the data, presented in Table 1, 
are estimated from the Monte Carlo simulation. Only statistical errors are quoted. In this data 
sample the only relevant contamination arises from r + r~ events with misidentified decay modes 
(r background). 

Table 1: Selection efficiencies and r background for the different decay channels, obtained from the Monte Carlo 
simulation and presented with statistical errors only. For this table all identified events are retained. 



r decay 


Efficiency (%) 


r Background (%) 


71-71 


57.57 ± 0.39 


24.18 ± 0.39 


71-7171° 


58.39 ± 0.19 


21.44 ± 0.18 


7T-7r27r° 


50.36 ± 0.31 


34.09 ± 0.34 


7T-37T 


54.29 ± 0.31 


16.42 ± 0.28 


7r7r -7T7T 


59.76 ± 0.19 


19.47 ± 0.17 


7171°-7t27T° 


52.12 ± 0.22 


31.92 ± 0.23 


7T7r°-37r 


54.66 ± 0.21 


13.96 ± 0.19 


7r27r°-7r27r° 


45.84 ± 0.50 


42.73 ± 0.56 


7r27r -37r 


46.98 ± 0.35 


27.69 ± 0.39 


37T-37T 


50.98 ± 0.48 


8.57 ± 0.36 



4.2 Tau direction of flight 

The reconstruction of the r flight direction is mandatory in this analysis to access the event 
observables, which are functions of the four-momenta of the final state particles in the r rest 
frame. This can be achieved in the semileptonic decays, for which the r direction lies on a cone 
around the total hadron momentum. For events with both taus decaying semileptonically the r 
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direction lies along one of the intersection lines of the two reconstructed cones in the case that the 
taus are produced back-to-back, with equal energies given by y/s/2, and m Vr = m 9T = 0. However, 
the two cones may not intersect due to detector effects or radiation. If the cones intersect, the 
event is considered twice using either solution. If the cones do not intersect, the particle momenta 
are fluctuated within their measurement errors and the event is accepted if the cones intersect in 
a minimum number of trials; the average direction is then used [9]. 

The effect of using both r directions has been studied with a Monte Carlo sample having 
approximately the same size as the data. Table 2 presents the statistical errors obtained in 
this analysis and using the correct r direction from the information at the generator level. Not 
distinguishing between the two r directions induces some degradation in the overall sensitivity for 
the four anomalous weak dipole moments. 

Table 2: Statistical errors obtained from a Monte Carlo sample approximately equal in size to the data sample, 
using this analysis and selecting the correct t direction from the information at the generator level. 





This analysis 


Correct r dir. 


<trc^ t [10- 3 ] 

^Im^ T [10~ 3 ] 


0.43 
0.76 


0.34 
0.58 


VRe dT [10- 3 ] 
Tlmd T [10- 3 ] 


0.39 
0.65 


0.36 
0.55 



4.3 Candidates and efficiency matrix 

The final selection for this analysis requires the r direction to be successfully reconstructed, as 
described in Section 4.2, and the event observables (W\, cos 9^, 0^, W2, cos^ 2 , <ph 2 ) to lie in 
their domains of validity. These requirements decrease the number of candidates by 21%, the 
main reason being the inability to reconstruct the r direction for some events. 

The final number of candidates in each decay topology is given in Table 3. The efficiency 
matrix e^- for i,j = it, p, tt2tt°, 3tt is calculated as a function of the generated polar angle cos 9^ 
separately in the barrel and endcaps; the dependence of on and is quite uniform and 
has been integrated out. Figure 3 shows the efficiencies for the barrel. The diagonal elements 
of this figure represent the identification of each decay mode, while the off-diagonal elements 
represent its misidentification. 

4.4 Detector effects 

The correct approach to introduce the detector effects in the likelihood formula would 
be via a smearing function depending on 12 variables: the set of event observables 
and the corresponding generated values. The indices i and j indicate the generated 
and the reconstructed channel, respectively. With the notation used here, Ty = 
Tij(Wi, cos6 hl ,(f) hl ,W2, cos6 h2 , (f) h2 , wf>, cos^f, 0^, cosflj^, 4>f ] 2 ). Because this function 
cannot be easily calculated, the detector effects are parametrised by the factorised smearing 
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Table 3: Number of final candidates in each decay topology and total number of events used in the analysis. 



Class 


Events 


Class 


Events 


71-71 


1901 


7r7r°-7r27r° 


6395 


7T-7T7T 


7844 


7T7T -7t37T 


5242 


7r-7r27r° 


2673 


7r27r°-7r27T u 


1125 


7T-37T 


2040 


7r27r u -37T 


1950 


7r7r°-7T7r 


8624 


3tt-3tt 


712 



Total number of events: 38506 



functions Dij(x, x^) (x = W, cos 6^, <f) h and i,j = rr, p, Tc2rr , 3tt). For x = W, cos6> fc , correlations 
are neglected, while for x = cos Oh, 4>h, the correlations are taken into account. 

The functions Dij(x, x^) give the probability that for generated i and reconstructed j the 
smearing introduced by the detector is (x — x^) for a certain generated x^ with reconstructed 
x. From the definition it follows that 

J Dij(x,x {0) )dx = 1 . (15) 

These functions are obtained with the SM Monte Carlo simulation by binning the (x,x^) plane. 
The binning has been chosen small enough to correctly convolve detector effects with the generated 
distribution. 

In the likelihood expression, Eq. 13, the detector effects are included by replacing the functions 
H? by 

H?(W,cos0 h ,<f> h ) = W ffr(wW,oosfli 0) ,^ 0) )£>y(W; W (0) ) (16) 

i J 

x A^cos^cos^A,^^^ 

The sum runs over all modes % which have been reconstructed as one of the modes j used in the 
analysis, whereby all possible r decay modes are included in i. The r branching fractions are 
taken into account implicitly in Eq. 16 because the full differential cross section (Eq. 12) contains 
the probability of generating a t + t~ pair decaying into specific decay modes with certain final 
state topologies. 

In terms of the effective functions iff the likelihood for each event reads 

Lij = E ^(^,^1^)^(^1,008^,^)^(^2,008^,^) . (17) 
H,u=0,...,3 

4.5 Calibration curves 

There are two sources for possible bias in the fitting procedure: (i) the detector effects are handled 
by the factorised smearing functions above which take correlations into account only partially, 
and (ii) radiative corrections are not included in the likelihood. It is thus necessary to evaluate 
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the adequacy of the fitting process. This is done with the SCOT Monte Carlo program [19], 
interfaced with TAUOLA for the r decays and with the full detector simulation. The program 
SCOT describes e + e~ — > t + t - production at an energy around the Z peak at tree level. It 
includes the anomalous weak dipole moments and all r spin effects. The initial state radiation is 
included by adding a simple radiator function [20]. 

The checks are performed by generating various Monte Carlo samples with different values of 
the anomalous weak dipole moments. The couplings a e , v e , a T , v T are set to their SM values. The 
anomalous weak dipole moments Re/i T , Im/z T , Re<i T , lmd T are varied one by one in an adequate 
region around zero. The dependence of the reconstructed values on the generated parameters 
is taken as linear. Significant deviations of the slopes from unity are found for certain decay 
topologies, and the offset for Re fi T is not consistent with zero for certain channels. These effects 
have been studied and are mostly related to not using the correct r direction and to background 
effects. In this analysis, this calibration for each anomalous weak dipole moment and decay 
topology is taken into account to obtain the corresponding individual measurements. The slopes, 
offsets and x 2 of t ne linear fits are presented in Table 4. 

The offsets and slopes were derived with a Monte Carlo program which includes only a first 
order radiator for the initial state bremsstrahlung. The offsets have then been verified with 
the KORALZ Monte Carlo program [21] interfaced with the full detector simulation. KORALZ 
describes e + e~ — > t + t~ production at an energy around the Z peak and with initial state 
bremsstrahlung corrections up to 0(a 2 ), final state 0(a) bremsstrahlung and 0(a) electroweak 
corrections. However, this program contains only longitudinal spin effects (i.e., the production 
terms used are R 00 , R 03 and R 33 ) and the anomalous weak dipole moments are set to zero. To 
check the effect of this approximation, a maximum likelihood fit was built to include the complete 
Roo, R03 and R33 terms and the anomalous parts of the other R^ v terms. The values obtained 
for the four anomalous weak dipole moments for each decay topology are consistent with the 
corresponding offsets computed with SCOT and the first order radiator. Table 5 shows the global 
differences between the offsets of the central values of the four anomalous weak dipole moments 
when using either KORALZ or SCOT. The final results of this analysis are obtained by correcting 
the individual measurements according to the KORALZ offsets and the SCOT slopes and by 
including the statistical error of the correction in the systematic uncertainty. The statistical 
error of the SCOT offsets is also included in the systematic uncertainty as these offsets have 
contributions from all R au terms. 



5 Systematic uncertainties 

The systematic uncertainty is calculated for each anomalous weak dipole moment and final state 
decay topology. The estimates are shown in Tables 6 to 9. The last row of these tables shows the 
combined systematic uncertainty from each source taking into account the correlations between 
channels. 

The ECAL effects are related to the uncertainty in the global energy scale and the nonlinearity 
of the response. The global scale is known at the level of 0.25%, through the calibration with 
Bhabha events. Global variations have been applied to each ECAL module and the effect is 
propagated to the fitted parameters. The nonlinearity of the response is related to the wire 
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Table 4: Results of the fit of the calibration curves for each of the decay topologies obtained with the SCOT 
program and a first order radiator, for [i T (top) and for d T (bottom). The slope a, the offset b and the \ 2 01 the 
linear fit are given for every case. The number of degrees of freedom is three for these linear fits. 



vjlldllllcl 


rve ^Li-r 


Im p T 




„ n n— 21 i, n n— 41 „ 2 

a [ID J o [ID J x 


a[io- 2 ; 


6[10- 4 ] x 2 


7T-7T 


7fi 4- c; 1/11_|_qo qoc 
(D./ it O.y 14.1 it o.o O.OO 


92.0 ± 6.3 -1.7 ± 4.3 0.39 


7T-D 

r 


93 3 ± 2 9 12 3 ± 1 7 2 28 


68.2 ± 3.6 0.1 ± 2.2 5.92 


7r-7r27T° 


104.7 ± 5.3 10.2 ± 3.2 0.09 


82.8 ± 6.7 0.5 ± 4.1 1.63 


7T-37T 


108.9 ± 5.4 5.4 ± 3.2 0.26 


77.3 ± 6.8 -0.4 ± 4.2 2.28 


o-o 


99.1 ± 3.1 9.3 ± 1.8 3.64 


60.4 ± 3.9 -1.1 ± 2.3 0.84 


p-7l27T° 


99.0 ± 3.9 0.1 ± 2.3 0.44 


69.8 ± 5.0 0.8 ± 3.0 4.45 


p-3iv 


89.7 ± 3.8 0.3 ± 2.2 3.55 


60.0 ± 4.9 0.9 ± 2.8 4.14 


■k2tt°-it2tc 


96.9 ± 9.9 -3.1 ± 5.7 1.72 


60 ± 14 -6.1 ± 7.9 0.28 


7r27T°-37r 


89.5 ± 6.6 -0.3 ± 3.9 1.75 


73.4 ± 8.9 -1.3 ± 5.3 1.10 


37T-37T 


102.0 ± 9.6 -9.6 ± 5.7 2.00 


55.0 ± 12.0 -4.0 ± 7.3 1.10 


i n annol 
vjlldllllcl 




Im d T 




^rin— 21 Ulin— 41 ^,2 

a [1U J o [iv \ x 


a[io- 2 ; 


b [io- 4 ] x 2 


7T-7T 


/O.Z it 4.1 — 1.0 it z.o O.zo 


51.3 ± 5.1 5.7 ±3.1 0.26 


n-p 


109.4 ± 2.9 1.2 ± 1.8 5.82 


90.5 ± 3.6 4.1 ± 2.3 4.31 


7r-7r27T° 


106.9 ± 5.4 1.3 ± 3.3 5.33 


91.4 ± 7.0 4.3 ± 4.4 1.17 


7T-37T 


113.5 ± 5.3 6.9 ± 3.3 1.27 


101.9 ± 7.0 2.2 ± 4.2 4.20 


p-p 


109.1 ± 3.2 -3.3 ± 1.9 1.83 


80.6 ± 4.1 0.6 ± 2.4 2.14 


p-7i2n° 


102.7 ±4.0 -1.1 ± 2.3 1.46 


70.4 ± 5.2 3.0 ± 3.1 3.08 


p-3lY 


96.4 ± 3.9 1.3 ± 2.3 3.76 


60.5 ± 5.0 3.3 ± 2.9 1.63 


7r27r u -7r27T° 


96.1 ± 9.9 0.5 ± 5.8 0.44 


75 ± 13 2.9 ± 7.8 1.08 


7r27r°-37r 


85.2 ± 6.8 4.3 ± 4.0 1.59 


66.4 ± 9.2 -1.1 ± 5.3 0.16 


37T-37T 


97.0 ± 10.0 2.0 ± 5.8 2.48 


60.0 ± 12.0 -2.1 ± 7.3 1.73 



Table 5: Global differences in the central values of the four anomalous weak dipole moments when using offsets 
from the KORALZ or SCOT programs. 



Parameter 


Shift 


Re p T [10" 3 ] 
Im p T [10" 3 ] 


-0.45 ± 0.21 
-0.43 ± 0.42 


Re d T [10" 3 ] 
Im d T [10- 3 ] 


0.05 ± 0.19 
0.81 ± 0.38 
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Table 6: Systematic uncertainties on Re \x T for the different channels. The last row gives the combined systematic 
uncertainty from each source taking into account the correlations between channels. The total systematic and 
statistical errors are shown in the last two columns. The values are expressed in units of [10 ]. The sources of 
uncertainty are explained in the text. 





ECAL 


TPC 


Align. 


rBF 


Wpar 


MC st. 


ai dyn. 


Fake 7 


0"sys 


Ostat 


7T-7T 


4.53 


0.84 


0.58 


1.55 


0.11 


14.30 


0. 


2.88 


15.39 


24.70 


TY-p 


2.90 


0.48 


0.04 


0.25 


0.12 


4.64 


0. 


0.23 


5.51 


9.35 


7T-7r27r° 


2.77 


2.98 


0.12 


1.12 


0.06 


7.52 


0.45 


0.54 


8.65 


14.81 


7T-37T 


0.63 


12.23 


6.95 


0.43 


0.14 


7.39 


1.23 


0.85 


15.98 


16.58 


P-P 


2.63 


1.23 


0.28 


0.37 


0.11 


4.43 


0. 


0.64 


5.36 


8.45 


p-n2n° 


3.20 


0.53 


0.42 


0.61 


0.06 


5.53 


1.88 


1.54 


6.90 


11.67 


p-3n 


1.55 


1.38 


1.94 


0.60 


0.08 


5.87 


0.61 


1.60 


6.77 


11.37 


7r27r°-7r27T 


10.24 


1.72 


0.86 


1.49 


0.14 


13.56 


11.37 


2.03 


20.69 


22.75 


7r27r°-37r 


4.50 


3.04 


2.38 


1.62 


0.03 


10.07 


2.91 


1.31 


12.22 


20.31 


37T-37T 


1.65 


7.82 


7.62 


0.89 


0.24 


13.17 


4.69 


4.61 


18.43 


25.55 


Combined 


0.72 


0.92 


0.43 


0.17 


0.08 


2.19 


0.43 


0.18 


2.60 


4.20 



Table 7: Systematic uncertainties on lm/j, T for the different channels. The last row gives the combined systematic 
uncertainty from each source taking into account the correlations between channels. The total systematic and 
statistical errors are shown in the last two columns. The values are expressed in units of [10 -4 ]. The sources of 
uncertainty are explained in the text. 





ECAL 


TPC 


Align. 


rBF 


Wpar 


MC st. 


ai dyn. 


Fake 7 


C"sys 


Cstat 


7T-7T 


1.60 


0.86 


0.42 


0.43 


0.09 


13.29 


0. 


0.97 


13.46 


25.06 


71- p 


0.86 


0.38 


1.52 


0.39 


0.05 


8.11 


0. 


2.44 


8.66 


16.60 


7r-7r27T° 


27.55 


3.40 


3.03 


1.58 


0.39 


14.92 


8.75 


3.24 


33.05 


29.65 


7T-37T 


11.23 


9.64 


9.98 


0.61 


0.27 


15.72 


8.94 


0. 


25.42 


32.96 


P-P 


3.29 


0.13 


0.85 


0.64 


0.09 


10.21 


0. 


1.09 


10.84 


18.09 


p-7r27r u 


32.05 


2.98 


1.13 


0.56 


0.06 


12.51 


2.01 


8.19 


35.57 


25.10 


p-3n 


9.60 


1.51 


0.57 


1.43 


0.08 


12.70 


2.04 


10.85 


19.40 


21.85 


TV ATI -7TA7T 


100.21 


2.20 


6.70 


5.49 


0.57 


35.31 


7.17 


10.17 


107.35 


62.38 


7r27T -37T 


17.50 


11.90 


5.75 


1.14 


0.22 


17.61 


12.21 


15.62 


34.43 


37.88 


37T-37T 


5.38 


5.70 


1.66 


0.99 


0.31 


34.50 


10.76 


17.05 


40.77 


64.46 


Combined 


4.27 


0.17 


0.30 


0.23 


0.05 


4.18 


1.12 


0.76 


6.10 


8.00 
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saturation constants of ECAL, which are different in the barrel and endcaps. These saturation 
constants have been fluctuated within their nominal errors while keeping the measured energy 
fixed at M z /2. 

The TPC systematic errors are related to the momentum measurement: (i) an effect due to 
the magnetic field acting similarly on positive and negative charged tracks, and (ii) a sagitta effect 
affecting oppositely positive and negative tracks. These two effects are calibrated with dimuon 
events and the corresponding corrections are applied to the r data. The systematic errors are 
then estimated by varying the corrections within their errors for each year. 

The systematic errors in the column labeled "Align." are due to a possible azimuthal tilt 
between the different parts of the detector. 

Variations in the r branching fractions are considered in the rBF column. This systematic 
uncertainty was determined from its effect on the calibration curves (Section 4.5). 

The experimental errors on the weak parameters sin 2 9w and Mz are propagated to the fitted 
values. Other weak parameters have negligible effect on the measurements. These effects are 
summarised in the column labeled "Wpar" . 

The finite Monte Carlo statistics also causes systematic uncertainties. The most relevant 
statistical uncertainty is for the KORALZ offsets. The statistical error of the offsets and slopes 
obtained with SCOT and the first order radiator (Table 4) are also taken into account. Finally, 
the statistical error in the calculation of the efficiency matrix is also considered. These effects are 
shown under the column "MC St.". 

Table 8: Systematic uncertainties on Red T for the different channels. The last row gives the combined systematic 
uncertainty from each source taking into account the correlations between channels. The total systematic and 
statistical errors are shown in the last two columns. The values are expressed in units of [10~ 4 ]. The sources of 
uncertainty are explained in the text. 





ECAL 


TPC 


Align. 


rBF 


Wpar 


MC st. 


a± dyn. 


Fake 7 


^sys 


0"stat 


71-71 


1.27 


0.58 


0.09 


0.77 


0.12 


8.33 


0. 


0.16 


8.48 


15.61 


Tl-p 


0.18 


0.40 


0.61 


0.35 


0.04 


3.91 


0. 


4.51 


6.03 


7.94 


7T-7r27r° 


8.38 


1.73 


1.58 


1.78 


0.09 


7.68 


3.36 


3.21 


12.63 


15.34 


7T-37T 


1.36 


0.57 


0.02 


0.27 


0.02 


7.22 


2.71 


1.00 


7.92 


15.05 


P-P 


0.53 


0.35 


0.49 


0.38 


0.04 


4.18 


0. 


0.17 


4.28 


8.04 


p-7r27r u 


1.45 


1.68 


0.63 


0.79 


0.08 


5.39 


0.79 


2.63 


6.52 


10.89 


p-3n 


5.37 


1.68 


0.79 


0.45 


0. 


5.56 


3.00 


2.26 


8.80 


10.89 


7r27r u -7r27r° 


8.75 


2.00 


1.10 


1.71 


0.07 


13.35 


8.18 


9.44 


20.47 


29.74 


7r27r u -37r 


5.26 


12.07 


10.96 


2.56 


0.05 


10.39 


4.88 


3.75 


21.11 


21.65 


3tt-3tt 


1.91 


4.75 


4.70 


1.64 


0.12 


14.06 


5.34 


2.25 


16.80 


29.39 


Combined 


0.41 


0.29 


0.26 


0.17 


0.02 


1.94 


0.90 


0.50 


2.30 


3.90 



The ai decay dynamics are not well described theoretically. The impact of this was evaluated 
in the past [22] by implementing several models in the analysis [23]. The implementation of those 
models is much more difficult in the present analysis. The uncertainty is estimated by means of 
three models: the Kiihn & Santamaria (KS) model [24] (used in the fitting formula), the Feindt 
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Table 9: Systematic uncertainties on Im<i r for the different channels. The last row gives the combined systematic 
uncertainty from each source taking into account the correlations between channels. The total systematic and 
statistical errors are shown in the last two columns. The values are expressed in units of [10~ 4 ]. The sources of 
uncertainty are explained in the text. 





ECAL 


TPC 


A 1 * 

Align. 


rBF 


Wpar 


MC st. 


a\ dyn. 


Fake 7 




Cstat 


7T-7T 


1.02 


0.59 


0.65 


1.25 


0.06 


14.95 


0. 


5.01 


15.87 


26.10 


ir-p 


o /to 
Z.4© 


U. /4 


n a 1 
U.41 


U.lo 


n no 


c nn 

o.yu 


U. 


U.y4 


^ .44 


13.00 


7T-7T Z7I 




fi ^fl 


9 7^ 


1 1 ZL 


U.ol 


1 ^ fin 


o.oy 


^ 9Q 
o.zy 


^4 7Q 


9Q ^ 
zy.oo 


7T-37T 


3.47 


3.38 


5.14 


0.79 


0.14 


11.90 


4.30 


8.03 


16.59 


23.50 


P-P 


2.91 


0.60 


0.55 


0.63 


0.04 


8.54 


0. 


0.26 


9.09 


15.07 


P-7V27T 


17.69 


0.44 


0.29 


1.20 


0.08 


12.46 


6.52 


2.54 


22.78 


24.64 


p-3n 


6.09 


9.99 


5.87 


0.36 


0.07 


12.40 


1.04 


5.91 


19.01 


23.31 


7r27r°-7r27r° 


9.86 


10.57 


17.41 


3.38 


1.15 


28.17 


40.10 


14.34 


55.96 


72.95 


7r27r u -37T 


8.12 


7.21 


3.72 


4.21 


0.11 


19.32 


1.73 


6.86 


23.93 


41.51 


37T-37T 


4.46 


12.78 


8.80 


0.67 


0.29 


30.10 


7.80 


4.19 


35.30 


59.22 


Combined 


5.74 


1.80 


0.70 


0.32 


0.04 


3.80 


0.79 


0.43 


7.20 


7.20 



model [25] and the Isgur, Morningstar and Reader (simplified) model [26]. The effects of the 
Feindt and IMR models on W and h 3 are calculated. The corresponding ratios with the W and 
the h 3 of the KS model are then used to scale the error. 

Another source of systematic error is due to fake photons generated by hadron interactions 
in the ECAL or by electromagnetic fluctuations. This quantity of photon candidates is 
underestimated in the Monte Carlo simulation compared to the data. This deficit was originally 
observed in a substantial disagreement between the data and the Monte Carlo simulation for the 
W distribution in the 7r27r° channel. This discrepancy has notably decreased after weighting the 
events of the Monte Carlo simulation according to the number of fake photons. This weighting 
was optimised for other r analyses [9]. Figure 1 compares the W distributions for the data and 
the Monte Carlo after the approximate weighting. In the end, the effect of fake photons is taken 
into account by removing fake photons in the simulation and using the difference in the fitted 
parameters as systematic uncertainty. 

The total systematic error and the statistical error for each channel and parameter are shown 
in the last two columns of the tables. 



6 Results and conclusions 

The final individual measurements of the four anomalous weak dipole moments are obtained 
applying the offsets and slopes described in Section 4.5. The results for the different decay 
topologies are presented in Figs. 4 and 5, including both the systematic and the statistical errors. 
All these measurements are consistent with the SM prediction. Figure 6 shows the relative weights 
of the different decay channels for the four measured anomalous weak dipole moments. 

The final combined results on the four anomalous weak dipole moments are listed in Table 10, 
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showing the statistical, systematic and total errors. The statistical correlations are given in 
Table 11. The final 95% CL upper limits derived from these measurements are presented in 
Table 12. 



Table 10: Final results on the real and imaginary terms of the anomalous weak dipolc moments. 



Parameter 


Fitted value 


Cstat 


&sys 


a 


Re/z T [10" 3 ] 
Im/i T [10~ 3 ] 


-0.33 
-0.99 


0.42 
0.80 


0.26 
0.61 


0.49 
1.01 


Red T [10" 3 ] ([HT 18 ecm]) 
lmd T [10" 3 ] ([10- 18 ecm]) 


-0.11 (-0.59) 
-0.08 (-0.45) 


0.39 (2.14) 
0.72 (4.00) 


0.23 (1.26) 
0.72 (4.01) 


0.45 (2.49) 
1.02 (5.67) 



Table 11: Statistical correlations between the fitted parameters. The individual correlations are presented in the 
off-diagonal elements. 





Re \i T 


Im fj, T 


Red T 


Im d T 


Re fx T 
Im \x T 


1.0 


0.006 
1.0 


0.028 
-0.055 


0.062 
0.034 


Red T 
Im d T 






1.0 


-0.003 
1.0 



Table 12: Upper limits derived from this measurement of the anomalous weak dipole moments (95% CL). 



Parameter 


Limit 




Re fj, T 
Im jj, T 


[io- 3 ] 
[io- 3 ] 


1.14 
2.65 




Re<i T 
Im d T \ 


[10- 3 ] ([10" 18 ecm]) 
[10" 3 ] ([10- 18 ecm]) 


0.91 (5.01) 
2.01 (11.15) 



These results supersede the previous ALEPH measurement of Red T [3]; the measurement of 
Re/x r , Im/i,,- and lmd T , presented in this paper gives the most stringent limits on these quantities 
to date. 
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Figure 1: The W observable for the p, tt2tt° and 3tt decays. The points are the data and the histogram is the 
simulation. Both distributions are normalised to unit area in each plot. Only the statistical errors are included. 
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Figure 2: The cos Qh observable for the four decay topologies. The points are the data and the histogram is the 
simulation. Both distributions are normalised to unit area in each plot. Only the statistical errors are included. 
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Figure 3: Efficiency function (cos ) in the barrel region. 
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Figure 4: Results on yu r for the various decay modes, including both the statistical and the systematic 
uncertainties. The results on Re[i T are shown at the top, and on Im^ T at the bottom. 
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Figure 5: Results on d T for the various decay modes, including both the statistical and the systematic 
uncertainties. The results on Red r are shown at the top, and on Imc? T at the bottom. The anomalous weak 
electric dipolc moment is assumed dimcnsionless in these figures. 
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Figure 6: Relative weights of the different decay topologies for the four measured anomalous weak dipole moments, 
normalised to the total weight. 



22 



